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Abstract. 

In this work we study the diffusion of non-interacting overdamped particles, 
moving on unbiased disordered correlated potentials, subjected to Gaussian white 
noise. We obtain an exact expression for the diffusion coefficient which allows us 
to prove that the unbiased diffusion of overdamped particles on a random polymer 
does not depend on the correlations of the disordered potentials. This universal 
behavior of the unbiased diffusivity is a direct consequence of the validity of the 
Einstein relation and the decay of correlations of the random polymer. We test the 
independence on correlations of the diffusion coefficient for correlated polymers 
produced by two different stochastic processes, a one-step Markov chain and the 
expansion-modification system. Within the accuracy of our simulations, we found 
that the numerically obtained diffusion coefficient for these systems agree with the 
analytically calculated ones, confirming our predictions. 
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1. Introduction 

The diffusion of particles on one-dimensional (ID) random media has been the subject 
of intense research. Since the introduction of the Sinai’s model it was recognized that 
disordered systems can exhibit normal and anomalous diffusion as well as normal 
and anomalous drift [1]. Several models of diffusion in disordered media have been 
introduced since then [2-4] in order to understand the underlying mechanisms leading 
to such behaviors. Particularly some simple deterministic models for transport in 
disordered potentials have been studied in order to understand the origin of the 
transport properties from the very deterministic dynamics [5-9]. A model for transport 
in random polymers that takes into account the pair correlations between monomers 
(to which we will refer to as the “particle-polymer model”) has been introduced by 
the authors in Ref. [9]. In such a work it was shown that the deterministic biased 
diffusion is influenced by the presence of the correlations in the polymer. Indeed, 
in that paper, the authors gave an exact expression for the particle current and the 
diffusion coefficient in terms of the first two moments of the “crossing times” and 
its corresponding pair correlation function [9], a result which is a consequence of the 
central limit theorem [10-12]. Based on the same arguments, it has been shown that 
when the particles along the polymer are placed at a finite temperature and subjected 
to a constant driving force, it is still possible to obtain an exact expression for the 
drift and diffusion coefficients in terms of the first and the second moments of the first 
passage time (FPT) as well as its pair correlation function [13]. In this work we will 
show that, for the particle-polymer model, the diffusion coefficient does not depend on 
the correlations between monomers along the polymer when the driving force is zero. 
The latter is a result that follows straightforwardly from the validity of the Einstein 
relation in disordered potentials at zero driving force. 

In consequence the paper is organized as follows. In Section 2 we introduce 
the model as well as some general results concerning the transport process previously 
established. Then in Section 3 we give the exact expression for the diffusion coefficient 
in disordered polymers for zero driving force and introduce a simple model for the 
potential to illustrate the independence of the diffusivity on the correlations. First we 
use a Markov chain to produce exponentially correlated polymers which are used to 
induce a correlated potential. Then, we perform Langevin dynamics simulations to 
estimate the diffusion coefficient of overdamped particles moving on such polymers. 
We also perform analogous numerical experiments in the case in which the polymers 
are produced by means of the expansion-modification system, which is known to have 
polynomially decaying correlations. Finally in Section 4 we give a brief discussion 
and the main conclusions of our work. Two appendices are given containing detailed 
calculations of our main results. 

2. The model 

We will consider an ensemble of Brownian particles with overdamped dynamics moving 
on a ID disordered potential V(x) subjected to an external force F. The equation of 
motion of any of these particles is given by the stochastic differential equation, 

^dXt = (/(W) + F)dt + godWt (1) 

where Xt stands for the position of the particle and Wt is a standard Wiener 
process. The constants F and 7 stand for the driving force and the friction coefficient 
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respectively. The constant po represents the square root of the noise intensity (or the 
thermal fluctuations), which according to the fluctuation-dissipation theorem satisfy 
the relation pg = Here f denotes, as usual, the inverse temperature times the 

Boltzmann constant, f = l/ksT. The function f{x) represents minus the gradient 
of the potential V{x) that the particle feels due to its interaction with the substrate 
where the motion occurs. As in Ref. [13], we model the substrate (or the polymer) as 
a chain of “unit cells” of constant length L. The unit cells represent the monomers 
comprising the polymer. We denote by A the set of possible monomer types, a set 
that is assumed to be finite or countable infinite. The polymer is represented by a 
bi-infinite symbolic sequence a := (... , a_i, oq, oi,...), where aj G A stands for the 
monomer type located at the j-th cell, for all j G Z. The set of possible random 
polymers will be denoted by JA according to the conventional notation in symbolic 
dynamics [14]. As in Refs. [9,13], we assume that the disordered potential V(x) is the 
result of the interaction of a particle with the random polymer. 

Let a; G K denote the particle position along the substrate a G A particle 
located at x feels a (random) potential that is not only function the particle position 
X, but also of the substrate, i.e. V(x) = ip(x,a). This is a consequence of the fact 
that the particle might interact with the whole polymer and not only with the closest 
monomer. If a; = 0 we assume that the particle is located at the beginning of the 
0-th monomer uq. In Ref. [13] it was shown that the random potential ip(x, a) has the 
following property, 


if(x + nL,a) = i/j[x,i7’^(a)j. (2) 

where cr : A^ —>■ A?‘ is the shift mapping defined as follows: if a, b G AA are such that 
(j(a) = b, then bi = a^+i for all f G Z. As in Refs. [9,13] we assume that the substrate 
is generated by some stochastic process through a stationary cr-invariant probability 
measure v on A^. This hypothesis is equivalent to say that the statistical properties 
of the polymer are translationally invariant. 

In Ref. [13] it has been shown that the particle flux Jeff and the diffusion coefhcient 
Heff can be written exactly as follows. 


Teff .— 
DeS ■ = 


lim 

t—¥(X) 


((M 

t 


L 

((t(O^L)))’ 


lim 

t—¥00 


Var(A,) 

2t 


2{{riO^LW 


( 3 ) 

( 4 ) 


where r(0 L) is a random variable, called the first passage time (FPT), defined 
as the time that the particle spends to reach, for the first time, the position Xt = L 
from the initial condition Xq = 0. The quantity gfj. is defined in terms of the variance 
of the FPT and its autocorrelation function as. 


C ? •— (('''^(0 ^ ^))) ~ (('''(0 —^ L ))) + 2 ^ Crim ). 

m—1 

where 

C'r(m) := (( t (0 — s- L)T{mL — >■ (m + l)j)))-((r(0^L)))^ 

It is important to stress that the notation ((•)) denotes a double average, one with 
respect to the noise and the other with respect to the disorder. In Ref. [13] it was 
pointed out that is important to distinguish between these two different averages. 
If we put an ensemble of non-interacting Brownian particles over a (fixed) random 
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polymer a we denote the average over the ensemble of particles as {■)n- This average 
will be referred to as the average with respect to the noise (w.r.n.). Once a certain 
observable has been averaged with respect to the noise, it still depends on the specific 
realization of the polymer a. Thus we need to perform a second average which should 
be carried out over an ensemble of different realizations of the random polymer. This 
average is taken with respect to a stationary measure v which characterizes the process 
producing the random polymers. This average will be referred to as the average with 
respect to the polymers (w.r.p.) and will be denoted by (•)p. Then, the double average 
{{•)) is actually achieved first by taking the average w.r.n. and after that by taking the 
average w.r.p., i.e., ((•)) = {{•)n)p- Additionally, we will use the following notation: 
Varn(Cl) := {0'^)n — {O)^ and Varp((!l) := {0‘^)p — {0)p denote the variance of the 
observable O w.r.n. and w.r.p. respectively. Along this line, Var(O) therefore denote 
the variance of the observable O with respect to both, the noise and the polymer 
ensemble, i.e., Var(O) := {{O'^)) — {{O))'^. 

Let us denote by ri(a) the mean FPT w.r.n., i.e., ri(a) := (t(0 L))^- Then 

the drift and diffusion coefficients can be written as [13], 

JeS = irp v ; (5) 

Lleff — T^noisy (6) 


where Unoisy and Udet are referred to as the noisy and deterministic parts of Heff 
respectively. These quantities are defined as follows 


Dr. 


Ddet — 


L2(Var„(T(0 ^ L)))^ 

’ 

L2Varp(ri) + 2L2 Crim) 

2(Ti)3 


(7) 

( 8 ) 


The autocorrelation of the FPT can also be rewritten in terms of Ti, which gives 


CA£) = (Ti(a)ri [a^(a)] )p - (Ti(a))J. (9) 

The reason for expressing the transport coefficients, Jeff and D^g, in terms of the 
mean and the variance (with respect to the noise) of the FPT is that such quantities 
can be written exactly in terms of quadratures. Indeed, in Ref. [13] it was proved that 
the mean FPT w.r.n., Ti, can be written as. 


Ti(a) = 7/3 ^ e '"^^^g+(a)g_[(T ""(a)] + 7 /?Jo(a), (10) 

m—l 


and the variance of the FPT, Varn[r(0 —>■ L)], as 

{ OO OO OO / 

EEE e-("+™+')^^^(?+(a)g+[a-"(a)] 

n—1 m—l 1 — 1 

X g_[a-”^-"(a)] 9 _[a-”-'(a)]) 
oo oo 

+ 2 EE e-(”+”^)^^'^g+(a)g_[a-"-'"(a)]/o[a-"(a)] 

n—1 m—l 
oo 

+ Ee-"^^^(7+(a)/i[a-"(a)] 

n—1 



Unbiased diffusion of Brownian particles on disordered correlated potentials 


5 


m—1 1 — 1 
oo 

+ 2 E e-^^^^q.[a-^{a)]h{a) + /4(a) 

m—1 

which defines completely the transport coefficients. Here, the functions q+,q- : 

M were defined as 

9 ±(a) := f dxexp(^±/3[tp{x,a) — xF]), (12) 

Jo 

and the functions B± : M x M and Q± : R x M were defined as, 



Q±{x,a) 



dyexp{± p[f{y,a) -yF]), 


B±{x,a) := exp ( ±/3[i/)(a;, a) — xF]). 


(13) 

(14) 


We also used the following short-hand notation for the involved integrals, 

/o(a) = / Q_{x,a)B+{x,a)dx, 

Jo 


A (a) 
/ 2 (a) 


[Q-{x, a)]^ B+{x, a)dx, 


Q+{x, a)B+{x, a)dx, 


r 

B+{x,a) / Q-{u,a)B+{u,a)dudx, 

Jo 

/ 4 (a) = f B^{x,a) f B+{u,a) [Q-{u,a)]^ dudx. (15) 

Jo Jo 

In Ref. [13] it was shown that, for the case of non-correlated random polymers, 
the transport of particles in this system arises an interesting phenomenon which is the 
non-monotonic behavior of the diffusivity on the temperature. In this work we will 
study the effective diffusion coefficient of overdamped particles moving on correlated 
and uncorrelated potentials. Bellow we will show with certain generality that the 
unbiased diffusion does not depend on the correlations of the polymer beyond the 
number of monomers that the particle interacts with. 



3. Unbiased diffusion on random polymers 

A remarkable property that can be inferred from our formalism, is that the diffusivity 
in unbiased disordered potentials does not depend on the correlations of the substrate. 
This means that the particles in disordered potentials are unable to feel the pair 
correlations between monomers comprising the polymer. Consequently, the transport 
properties of a system of particles diffusing on correlated polymers are equal to those of 
a system of particles moving on an “equivalent” ensemble of uncorrelated polymers. 
Recently [15] this kind of behavior has been shown to occur in a similar model of 
diffusion on random potentials. In Ref. [15] Goychuk and Kharchenko calculated 
the diffusion coefficient of overdamped particles diffusing on a correlated random 
potential with a Gaussian distribution. Surprisingly their result does not depend on 
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the correlations, a conclusion that is reached by using the Einstein relation. Here we 
prove that for our model of random potentials the diffusion coefficient is independent 
on the correlations if the substrate is ergodic and has decay of correlations. The word 
ergodic in this context means that the polymer average can be taken as a time average 
if we think of the polymer as a discrete-time stochastic process. 

Bouchaud and Georges in Ref. [4] have shown that the Einstein relation, /r = 
fiDes{F = 0), remains valid for the class of disordered potentials we are considering. 
The constant p. stands for the mobility at zero bias, p = limF-s-o From the 

Einstein relation it is easy to obtain the diffusion coefficient since 

Des{F = 0) = - lim ^ . 

p F^o F 

In Appendix A we prove that the first moment of the FPT the divergence occurs as 
F~^. This implies that the particle current vanishes, as expected, in the limit F ^ 0. 
Actually the particle current behaves linearly in a small neighborhood of E = 0. This 
allows us to obtain the effective diffusion coefficient as follows 


DefiiF = 0 ) 


iPZ+ZJ 


(16) 


Here the constants Z± are defined as 


:= {g±(a))p = 


exp [±/3'i/)(a;, a)] dx 


Notice that the effective diffusion coefficient depends only on the quantities Z± 
which are local averages^ in the sense that the polymers involved in the averages are 
not shifted with respect to each other (as it occurs in the correlation function given in 
Eq. (9)). Thus, the expression (16) for the effective diffusion coefficient implies that, 
on ergodic polymers with decay of correlations, the diffusivity does not depend on the 
correlations between monomers at distant sites if the particle only interacts with the 
closest monomer to it (i.e., if the potential 'tp{x,a) only depends on one “coordinate” 
of a). This means that particles on polymers with long range correlations exhibit 
the same diffusivity as the particles that moves on polymers without correlations but 
with the same one-dimensional marginals. More generally, the diffusivity of a particle 
moving on a random polymer will be influenced by the pair correlations between 
monomers as distant as the range of the potential induced by the interaction. Here, 
the range of a potential stands for the maximal distance between two monomers having 
a non vanishing interaction with the particle. This makes the averages with respect 
to the polymer ensemble to depend on the range of the interaction. In other words, 
if the range of the potential 'f{x, a) is n, then, ip can be seen as a function of n 
coordinates of a, i.e., ipix, a) = ipiy, aj, aj+i ,..., aj+n-i)- Therefore, to perform the 
polymer averages Z± it is necessary to have at least the n-dimensional marginals of 
the probability measure v defining the polymer ensemble i.e., we need to know the 
joint probability distribution of n monomers, P(aj, aj+i ,..., aj+n-i)- With exception 
of the Bernoulli measure, the n-marginals cannot be decomposed, in general, as the 
product of the 1-dimensional marginals. In the following we give some examples to 
better illustrate these situations. 
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In order to exemplify the independence on the correlations of the diffusion coefficient 
at zero tilt strength we use the potential model introduced in Ref. [13]. The potential 
model V{x) used in such a work is of range 1 since the potential only depends on 
the monomer at which the particle is located. Explicitly the potential is defined as 
follows. If X represents the particle position and we write as x = y + nL, then 


V{x) = -fix, On) 


any if 0 < y < L/2 
an{L-y) if L/2<y<L. 


(18) 


Here n € Z stands for the index of the unit cell at which the particle is located and y 
is the relative position of the particle along such unit cell. The varible a„ is a random 
variable representing the monomer type at the nth unit cell. This potential model is 
shown schematically in Fig. 1. We observe that such potential is symmetric on every 
unit cell with a maximum located at y = 1/2. The height of the potential assumed to 
be random taking values from a finite set. Since the height of the potential is given by 
anT/2 we can assume that a„ represents the random variable, for every n G Z, which 
can take values from a set A := {fj : I < j < k}. 



Figure 1. Schematic representation of the potential model, (a) The potential 
profile on the 0th unit cell, (b) A realization of the random potential with a 
few unit cells. In this case every monomer along the chain can be taken among 
three possible monomers {k = 3) with heights Vi = fiLl2, V 2 = / 2-^/2 and 
V 3 = fsLl2. To perform analytical calculations as well as numerical simulations 
we have taken the values fi = 0.8, /2 = 4.2, and /$ = 9.0 (see text). 

We should stress here that the sequence a :=(..., a_i, oq, oi, ■. ■) G A?' represents 
the polymer (with a„ the corresponding monomers for n G Z). The values fj represent 
the “slopes” that the potential can take and, in some way, stand for the possible 
monomer types from which the polymer is built up. It is clear that the proposed 
potential depends only on one monomer, i.e., V{x) = ipiy, an) it x = y + nL. Since we 
are considering the Bernoulli measure on A^, we only need to specify the probability 
that a given monomer a„ equals a monomer type fj for 1 < j < k. In other words, 
we need to specify the I-dimensional marginal, i.e., jz(a„ = fj) =: p(fj) to perform 
the polymer average. With these quantities we can state explicitly how to take the 
average with respect to the polymer ensemble. If ft, : ^ > K., then we have 

k 

(ft(a))p = '^h{fj)p{fj). 

j=i 


( 19 ) 
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Notice that this average does not depend on n, which reflects the fact that the chosen 
measure is translationally invariant (or shift-invariant). With this potential model we 
have that the integrals Z_|_ and Z_ can be done exactly, giving 


= 2 


exp (±/3aL/2) — 1 
±f3a 




i=i 


±/ 3 /. 


( 20 ) 


The presence of correlations does not change the above average. This is because the 
function to be averaged only depends on one coordinate of a G Thus, to perform 
the average over the polymer ensemble we only require the 1 -dimensional marginals. 
Thus if two stochastic processes producing the (random) substrate have identical 
1 -dimensional marginals then they will exhibit the same diffusion coefficient at all 
temperatures. To confirm that it is the case, we simulate the Langevin dynamics of 
10 000 overdamped particles diffusing on several random polymers. Then we calculate 
the diffusion coefficient from the exact expression given in Eq. (16) and compare it 
with the diffusivity obtained from the simulations. To perform the referred numerical 
simulations first we built the polymers using a set of fc = 3 monomer types. The 
slopes defining the corresponding monomer types are chosen as /i = 0.8, /2 = 4.2 and 
/a = 9. The set of possible monomer types will be represented by A := {1, 2, 3}. Thus, 
a realization of the polymer a = (... ,a_i,ao,ai,...) has an associated sequence of 
potential profiles with “slopes” (..., /o_n faoj fan ■ ■ ■)• In. our numerical simulations 
we only vary the way in which the polymer is built up. For numerical reasons we fix 
the parameters L and 7 to one throughout the rest of this work. 

In Fig. 2 we observe the diffusion coefficient, as a function of the temperature, 
for particles moving on random polymers built up by means of a Bernoulli process 
and a Markov chain. On one hand, the Bernoulli process is defined uniquely by the 
1-dimensional marginals, pj = v(an = fj) for j = 1,2,3. Notice that the latter does 
no depend on the lattice site n, which means that the measure is shift invariant. For 
our numerical simulations we use the equidistribution for the monomer types, i.e., 
Pi = P 2 = P 3 = 1/3. This means that the polymer is built up at random as follows: 
at each unit cell we chose at random a monomer type from the set A := {1, 2, 3} with 
the same probability. Such a choice is independent of the monomers chosen in the 
rest of the chain. On the other hand, the Markov process that we use to built up the 
polymers is defined through the stochastic matrix P : ^ x ^ [0,1] given by 

/ 0 p q \ 

P= g 0 p , (21) 

\ p q 0 J 

where p and q are parameters such that q = 1 — p. It is easy to see that this matrix 
is doubly stochastic and the unique invariant probability vector tt = irP is given 
by TT = (4, 4, 4). It is known that the correlations for any Markov process decay 
exponentially fast [16]. Contrary to the Bernoulli process, the monomer type that is 
chosen at random at a given cell depends on the choice already made at an adjacent 
unit cell. Notice that although the Bernoulli and the Markov processes are different 
in nature as explained above, they have the same 1-dimensional marginals. This 
means that the probability that a given monomer type appears along the polymer 
is the same regardless the process (Markovian or Bernoulli). This is a consequence 
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Figure 2. Effective diffusion coefficient for overdamped particles on “Markov” 
and “Bernoulli” polymers. We simulate the dynamics of 10000 overdamped 
particles on random polymers built up by means of a Markov chain. The stochastic 
matrix corresponding to such a Markov process is given in Eq. (21) having a 
parameter p. We calculated the diffusion coefficient from the simulations for three 
different cases p = 0.50, p = 0.70 and p = 0.90 of the Markov chain. We also 
calculate the diffusion coefficient for the case of the Bernoulli process. We observe 
that the diffusion coefficients for the Markov chain do not vary with the parameter 
p and coincide with the diffusion coefficient obtained from the Bernoulli process. 
Moreover, these diffusion coefficients are consistent with the values predicted by 
Eq. (16). This gives an example of the fact that the diffusion coefficient (in absence 
of driving force) is independent on the correlations present in the polymers. 


of the fact the components of the invariant probability vector tt = (1/3,1/3,1/3) 
of the stochastic matrix P, actually correspond to the 1-dimensional marginals of 
the stationary Markov chain. In Fig. 2 we can appreciate that, within the accuracy 
of our numerical simulations, the diffusion coefficient for different realizations of the 
polymers coincide, regardless if they are built up by means of the stationary Markov 
chain or by means of the Bernoulli process. Moreover, these numerical experiments 
fits satisfactorily (within the accuracy of our simulations) to the exact curve calculated 
from Eq. (16). 

3.2. Potentials of range n = 2 

In order to test the independence on long-range correlations of the diffusion coefficient 
we simulate the dynamics of overdamped particles on random polymers built up by 
means of the expansion modification (E-M) process. The E-M system is a stochastic 
process that was introduced as a simple model exhibiting spatial 1// noise [17,18] and, 
in particular, it was used to understand the long-range correlation present in DNA 
sequences [19,20]. Such a system is defined by means of two fundamental processes 
as follows. Consider a “seed” (a symbol in the binary alphabet {0,1}). The seed is 
expanded with a probability p and it is modified (it is replaced by the complementary 
symbol) with probability 1 — p. In other words, if a € (0,1} then, under the E-M 
process this symbol is subjected to the transformation, 

J a with probability 1 — p 

1 aa with probability p. 
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and these transformations are extended coordinate-wise to words of arbitrary length 
in order to define the infinitely iterated E-M process. Here a stands for the 
complementary symbol, i.e., 1 = 0 and 0=1. This process allows the seed to grow in 
such a way that the resulting string of symbols becomes infinite with probability one. 
This stochastic process has been studied within the context of symbolic dynamics in 
Ref. [21]. In that work it was rigorously proved that the E-M system has a unique 
stationary measure having polynomial decay of correlations for an open set of values of 
the parameter p. Specifically, they proved that the pair correlation function C(i) has 
a behavior of the form C{£) x where 9 is an exponent depending on the parameter 
p, i.e. 9 = 9{p). The authors of Ref. [21] also conjectured that this behavior for the 
correlation functions occurs for almost all the values of the expansion probability and 
gave an explicit formula for the corresponding exponent 0(p), 


9{p) = 


log(l +p) -log(|2p- 1|) -log(|3p- 1|) 


( 22 ) 


log(l-kp) 

The treatment of the E-M system by means of symbolic dynamics allows us to obtain 
exact formulas for the n-dimensional marginals defining the (cr-invariant) stationary 
probability measure [21]. Indeed, the n-dimensional marginals are obtained as the 
stationary probability vectors of a hierarchy of stochastic matrices of size 2" x 2". 
Every stochastic matrix in such a hierarchy gives the “transition probability” from 
one word of size n to another one under the action of expansions and modifications 
on the polymer. Following the treatment of Ref. [21], we prove in Appendix B that 
the stochastic matrices for the cases n = 1 and n = 2 are given by. 


Ml = 


M2 = 


1 -l-p 


l+p 


2p 

1-p 

2p^ 


1-p 

2p 


V 


+ pq 

p^ + 2pq 

pq 


pq 

p 2 

q^ 


pq 

q^ 

p 2 


q^ 

pq 


pq pq 


P - 

2 p^ 


The respective invariant probability vectors vi 
given by 

'1 V 
2 ’ 2 , 


A 


2pq 
pq 

[0,1] and V 2 


(23) 


(24) 


+ pq ) 


■^1 = 7:: 7: 


[0,1] are 


(25) 


V2 = 


3 — 2p 1—p 1—P 3 — 2p 


10 - 8p’ 5 - 4p’ 5 - 4p’ 10 - Spy ’ 


(26) 


which are obtained as the left eigenvectors associated to the largest eigenvalues of Mi 
and M 2 respectively. With these quantities we can evaluate the average w.r.p. of any 
function on that depends on one or two coordinates of the polymer. Specifically, 
if we have a function depending on one coordinate of the polymer, p : A ^ K., then 
the average w.r.p. is defined by 

1 

{9{o))p = '^9{a)vi{a). (27) 

a=0 


On the other hand if we have a function depending on two coordinates of the polymer, 
—>• R, then the average w.r.p. should be achieved as follows, 

1 1 

{h{a,b))p = EE h{a, b)v 2 {a, b), 

a=0 b—0 


( 28 ) 
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Figure 3. Effective diffusion coefficient of overdamped particles moving on E- 
M polymers with expansion probability p = 0.80. In this figure we display the 
effective diffusion coefficient estimated from the Langevin dynamic simulations 
(filled circles) and the theoretical prediction from Eq. (16) (solid line). This figure 
clearly shows that the effective diffusion coefficient obtained from the simulations 
is consistent with the theoretical prediction, at least within the accuracy of our 
numerical experiments. 


a summation that will be represented alternatively as 

{h{a))p= ^ h{a)v2{a). (29) 


Now let us define a toy model for the potential induced by the particle-polymer 
interaction. We will assume that the particle interacts with two monomers when it 
is placed on the polymer. This interaction potential will depend on two coordinates 
of the random polymer. We will assume additionally that the induced potential i/' is 
piece-wise linear. Given a particle placed at the nth unit cell on the polymer a, we 
write the particle position x as x = y + nL with y G [0,1). Remember that in this 
case the polymer is built up from a binary “alphabet” A = {0,1} (since we have only 
two symbols in the expansion-modihcation process). Then the potential is defined as 


tp{x,an,an+i) 


ip{an,an+i)y if 0<y<L/2 
(p{an,an+i){L -y) if L/2<y<L. 


(30) 


The function ip : A? y. A? ^'M. appearing in the equation above is such that, 

( fi if a = 0 and 5 = 0 

ip{a, 5) = < /2 if (a = 0 and 5=1) or (a = 1 and 5 = 0) (31) 

( /a if a = 1 and 5=1. 

The values fi for i = 1, 2,3 are chosen as in the preceding example, /i = 0.80, /2 = 4.2 
and /a = 9.0. With this model we have that the resulting potential on a given unit cell 
depends not only on the monomer type located in such a unit cell but also depends 
on the right adjacent monomer type. Note that the potential profile defined by (30) 
is the analogous to the potential profile defined by Eq. (18). The difference is that in 
the present model the slopes are defined by two coordinates of the polymer, instead of 
one as in the potential model of Eq. (18). With these definition we can now evaluate 
the quantities and Z_ defined in Eq. (17) which allows us to obtain an expression 
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Figure 4. Effective diffusion coefficient of overdamped particles moving on E- 
M polymers with expansion probability p = 0.90. In this figure we display the 
effective diffusion coefficient estimated from the Langevin dynamic simulations 
(filled circles) and the theoretical prediction from Eq. (16) (solid line). We should 
notice a clear (non-random) deviation of the simulations with respect to the 
theoretical curve. In this figure we also display the effective diffusion coefficient 
of particle on the polymers produced by the E-M process but after shuffling the 
symbols. The shuffling process “destroy” the correlations and, as we can see, this 
diffusion coefficient differs from both, the theoretical result and the diffusivity 
originally calculated in the correlated polymer. This result evidence the fact 
that the E-M process has not reached its stationary state and the fact that the 
diffusion process in long-range correlated polymers undergoes a long transient 
resulting possible in a very slow convergence of the diffusion coefficient. 


for the diffusion coefficient through Eq. (16). Notice that the integral appearing in 
Eq. (17) is the same as in the case of potentials of range n = 1. Thus, it is easy to see 
that, 


= 2 


exp (±/3(^(a)L/2) — 1 


= 2E 

slGA'^ 


±/3(f{a.) 
exp {±l3ip{a.)L/2) — 1 


±/3(/?(a) 


U2(a). 


(32) 


In Fig. 3 we show the effective diffusion coefficient for particles along an E-M 
polymer with expansion probability p = 0.80. In this case the pair correlation function 
of the polymer is approximately given by C(i) w £^( 0 -S 0 )^ where 0(0.80) Rs 1.296, which 
implies that the correlations decay too slowly. To perform the numerical experiment 
(from which we estimate the effective diffusion coefficient) we simulated the Langevin 
dynamics of 35 000 particles, every particle moving on different realizations of the 
polymer, during a total time of 10® arb. units. We can appreciate in Fig. 3 that even in 
this situation, where the polymer has long range correlations, the diffusion coefficient 
obtained from the Langevin dynamics (fifled circfes) fits satisfactorily (within the 
accuracy of our simulations) with the diffusion coefficient predicted by Eq. (16) (solid 
line). This means, again, that the diffusivity does not depend on the large correlation 
present in the disordered substrate. 

We now test our formula (16) for the effective diffusivity in a more dramatic 
situation. This is the case where the correiations decay much more slowfy in 
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such a way that the correlation function behaves asymptotically as C{£) « £^(0-90) 
where 0(0.90) w 0.521. In this case, a typical realization of the E-M polymer has 
large ordered domains which are responsible of the large correlations present in the 
(disordered) substrate. It is natural to expect, in the present case, large deviated 
results from the theoretical ones. This is, effectively, the case as we observe in 
Fig. 4. The effective diffusion coefficient obtained from Langevin dynamic simulations 
(filled circles) deviates moderately (but with a clear non random tendency) from the 
theoretical curve predicted by Eq. (16) (solid line). This behavior could be the result of 
two phenomena naturally present in the system. First we should remind that the E-M 
process relaxes slowly to its stationary state, especially for large values of the expansion 
probability p [21]. As it was stated in Ref. [21], we should be aware of this phenomenon 
since it could lead to distorted observations. On the other hand, another phenomenon 
that could also be responsible of the deviations observed in Fig. 4 is the fact that the 
diffusion of the particles is by itself a very slow process. Indeed, the effective diffusivity 
is a property of the whole polymer, and if the particles are unable to visit the whole 
polymer, the diffusion coefficient would have large statistical errors. To verify that 
these two phenomena are actually present in the numerical simulations we performed 
the following experiment. We toke the polymer produced by the E-M process (with 
p = 0.90) and we shuffled its monomers in order to “destroy” the large correlations 
on the substrate. This procedure was done by taking two monomers at random and 
interchanging its position along the polymer. We performed 10® “interchanges” on 
every polymer of 10® symbols long. Then we placed the particles on the shuffled 
polymers and performed the corresponding Langevin simulations. In Fig. 4 we display 
the diffusion coefficient calculated in this way (open squares). We appreciate that the 
resulting diffusion coefficient still deviates from the theoretical one. This means that 
the deviation of the diffusivity observed in the long range correlated polymers does 
not come from the correlations at all, but from the fact that the polymer itself, has not 
reached its stationary state. This phenomenon can be appreciated in Fig. 5. In such a 
figure we plot the relative frequency of occurrence of the word 00 for 50 realizations of 
the random polymer of 10® symbols long. We can observe that for the case p = 0.80 
such relative frequency (filled circles) coincide in the average with the theoretical value 
(solid line), the two-dimensional marginal V2(00), predicted by Eq. (26). On the other 
hand, for the case p = 0.90, we can appreciate that the relative frequencies obtained 
from the simulations (stars) deviate largely from the theoretical prediction (dashed 
line) compared to the case p = 0.80. Indeed, the standard deviation of the relative 
frequency for the case p = 0.80 is An w 0.001550471 while for the case p = 0.90 
we obtain An « 0.02874428 which differs in one order of magnitude. This behavior 
clearly affects the effective diffusion coefficient obtained from the simulations even in 
the case where the correlations are destroyed by shuffling the symbols. 

4. Discussion and conclusions 

In this work we have considered a model for unbiased diffusion of particles moving 
along disordered correlated polymers. We have proved that if the particle only 
interacts with the monomer closest to it, the diffusivity does not depend on the 
correlations between monomers along the polymer. In general we have proved that 
if the particle interacts with a few monomers on the polymer, it will depend only on 
the correlations of monomers as distant as the range of the interaction. This result 
was proved under the assumptions (z) that the Einstein relation is valid in this kind of 
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Figure 5. Two-dimensional marginal of the polymer corresponding to the word 
00. We display the relative frequency of occurrence of the word 00 along every 
realization of the polymer. Every realization has 10® symbols and we produced 
50 different realizations for the cases p = 0.80 and p = 0.90. We observe that in 
the case p = 0.80 the numerically obtained two-dimensional marginal U 2 ( 00 ) is 
consistent with the predicted by Eq. (26). We also appreciate that for the case 
p = 0.90 we have large statistical deviation from the theoretical value compared 
to the case p = 0.90. Actually, the standard deviation for p = 0.90 is an order of 
magnitude larger than the standard deviation of the case p = 0.80. 


systems and (ii) that the polymer, regarded as a discrete-time stochastic process, is 
ergodic with decay of correlations. We illustrated the independence of the diffusivity 
on the polymer correlations by means of a simple model in two different cases. In the 
first case we produced the polymers by means of a three-states Markov process. Then 
these polymers were used to build up the potential felt by the particles. In this case 
we observed that the correlations does not alter the diffusivity and the corresponding 
effective diffusion coefficient that we obtained from the simulations was consistent with 
the exact result we derived. In the second case that we analyzed, the polymers were 
produced by means of the E-M process. The E-M process is known to exhibit long 
range correlations and the corresponding scaling laws for the correlation functions are 
known exactly. This feature of the E-M process allowed us to study the diffusivity in 
a system with a very slow decay of correlations. We could verify the independence 
on the correlations of the diffusivity; however we found that the effective diffusion 
coefficient seems to deviated moderately in the case in which the correlations decay 
as C{£) « £S(o.90) 0(0.90) « 0.521. We argued that this apparent inconsistence 

might be due to the fact that the E-M process relaxes slowly to its stationary state. 
Moreover, we observed that the numerically obtained two-dimensional marginal W2(00) 
(corresponding to the realizations of the polymers used for the numerical simulations) 
exhibits large statistical deviations on every realization of the polymer. This fact 
clearly affects the diffusivity along every polymer, resulting naturally in a deviation 
from the theoretical result. However it is also clear that the diffusion process by itself 
could be a very slow process since the particles should explore all the polymer in 
order to have an accurate estimation of the diffusivity. However, due to the existence 
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of large ordered domains in the polymer we could face a very slow convergence of 
the effective diffusion coefficient. It is therefore necessary a more detailed study of 
this class of systems particularly focusing in understanding how long is the transient 
and how does it behave. Knowing these properties could shed some light about the 
diffusion processes occurring in real systems in which the long-range correlations are 
naturally present, such as the DNA. 
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Appendix A. Diffusion coefficient in the limit of zero tilt strength 


In Ref. [13] it has been shown that in the polymer-particle interaction model for 
diffusion on random media the particle flux can be written exactly in terms of the 
mean FPT, 


'^efT — 


{Tih 


where the mean FPT w.r.n. Ti is given by 

OO 

Ti(a) =^I3Y^ e“™'’^^g+(a)g_[cr"’”(a)] 4-7/3/o(a). 


(A.l) 


(A.2) 


m—1 


Our goal here is to prove that in the limit of A = 0 the mobility /i, 

, J{F) 

is well defined for any stationary measure ly defining the polymer ensemble. The only 
condition we impose is the decay of correlations even at a very slow rate. In other 
words, we impose a mild condition on the rate of mixing of the dynamical system, 
which in turn implies ergodicity [22]. By ergodicity we mean that the time-average 
of any regular observable O : —>■ R, can be interchanged by its space average with 

respect to the shift-invariant probability measure ly, i.e., 

1 n-l 

/ O dv = lim — O o aUa) almost everywhere. 

Ja^ n^oo n ^ 

Here our observables will be q± and Ti which are non-negative and for any potential 
■0 bounded by above, they are also bounded, and furthermore are i.e. 


{Ti)p = [ Tidiy = f \Ti\diy<oo. 
Ja^ Ja^ 


Next, we will obtain an expression for the asymptotic behavior of Jes/F as F 
approaches to zero. In order to do this, let us first estimate bounds for (ri)p. We 
assume that the system {A^,a,iy) satisfies that for every / G L°°(^) and g G U^(iy) 
there exists a decreasing-to-zero function r(m) such that 


f ■ go a”^dv - / fdv / gdv 


< ll/lloo||g||i •r(m). 
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where || • ||oo and || • ||i denote the and norms respectively. The rate 

function r(m) will be precisely given below. Let us write C := ||/||oo||3l|i) this is a 
constant since q± are essentially bounded and L^(v), and thus we are able to give the 
following bounds 


(^i)p ^7/3 X! ® (^y J ± Cr{ni)j j 


^ 7/3 


(g+)p(g-)p 
1 - e-dPL 


± C £ + f hdi^] . 

m^l ) 


(A.3) 

(A.4) 


Which makes that the particle flux (Eq. A.l) be bounded by above and below as 
follows, 

L 

< _ 

7/3 T C Em=i + / hdv^ 

Next, let us consider the approximation as F —s- 0 of the following quantity. 


lim F 

F^O 


( (g+)p(g-)p 
^ 1_ e-PPL 


m—1 



The third part goes to zero since / I^dv is constant. By limit identities, the first part 
is given by 

lim ^ ^ = AA. 

F^o 1 - e~d^^ jdL 

It remains to estimate the summation involving the rate function, for which we make 
use the assumption on the decay of the correlations. Let us assume that the rate 
function r{m) = m~°‘ for any a > 0. This implies that the sum 


OO 

m—1 


can be written as the polylogarithm function [23], defined as, 

OO ^ 

Li„(x):=^|^. (A.5) 

k^l 

It is known that the polylogarithm fnnction has the asymptotic behavior [23] 

Lh(e-«) « r(l - s)r-i (A.6) 

for 1^1 <C 1 and for non-integer s (in particular for 0 < s < 1), being T the well known 
gamma function. If we consider a Fq > 0 small enough such that PFqL I, for any 
fixed [3 and L, then, one may write for every F < Fq, 

OO 

FC e-^dFL^-a _ (jp _ _ a){l3FL)°‘-^, 

m—1 

for all a > 0. In this way we have that 

OO 

lim FC = lim C • r(l - a)(/3F)“-^F“ ^ 0, 

F^O ^ F^O ^ 



Unbiased diffusion of Brownian particles on disordered correlated potentials 


17 


for all a > 0 and I3,L > 0. Now, using the Einstein relation we see that the effective 
diffusion coefficient is given by 


DeB = lim 


1 

7/3 ’ 


which is the expression for the diffusivity anticipated in Eq. (16). 


On one hand, the above treatment gave us the desired expression for the diffusivity 
when one takes firstly the average on the state space and after that the approximation 
when F goes to zero. On the other hand, one would expect that the expression (16) 
remains valid when one takes first the approximation E ^ 0 and after that the average 
with respect to the polymer ensemble. Here we give a result that enables us to 
argue that this could be actually true. Nevertheless, we left for further research the 
mathematically rigorous treatment. The importance on this kind of result is that the 
Einstein relation should hold even if we do not perform the average over the polymer 
ensemble. This means in particular that the mean FPT Ti should behave as F~^ for 
E ^ 0 which in turn implies Eq. (16) as we will see below. 

Eor the sake of convenience, we introduce a quantity of interest. Let N G N and 
define the observable gN ■ —>■ R as, 

- OO 

5 ^(a) := ^ E o a™(a), (A.7) 

m—l 

for some z G M“''. The discussion in the rest of this appendix is to justify the 
approximation 

Ee—^/^/oa-(a)«- / / = ^(/(a))p, 

for large N. 

For the moment, let us go to expression for the mean FPT and see that the 
approximation above may be applied to the summations appearing in Ti in Eq. (A.2) 
if we identify N with 1/F, 


Ti(a) = 7/3?+(a) E e ""(a)] + 7/3/o(a), 

m—l 

« 7/3g+(a)-g^(g-(a))p +7/3/o(a). 


(A.8) 


Notice that in the limit E —> 0 the mean FPT goes to infinity as F~^. The term 
7/3 Jo (a) appearing in the above expression remains finite if E = 0, and therefore it 
can be neglected in further calculations. Then we have that the particle flux behaves 
asymptotically as 

L FL^ 1 

"E (<?+(a))p(<7_(a))p’ 

when F is small enough. Now, using the Einstein relation we see that the effective 
diffusion coefficient is given by 

1 


D^e = i lim 


P F^o F -fPZ+ZA 
where are given by Eq. (17) and which is the expression for the diffusivity 
anticipated in Eq. (16). 

Next, in order to justify our claim, we give the following proposition. 
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Proposition 1 Let {A^, a, B, v) he a measure-preserving dynamical system, where 
a : A^ is the shift mapping, B is a Borel sigma-algehra on A^ and u is a 

shift-invariant probability measure satisfying decay of correlations at rate for 

some a > 0. Let f : A^ —5- R 6e o square integrable function on the full shift , i.e.. 


lA^ 


f^dv < oo. 


And let gN he the observable defined in (A.7). Then, for every e > 0 there exists a 
sufficiently large N such that 


a G 


gN{s) [ f dv 
^ d A^ 


> e 


< 


C-r{l-a)z° 

lV“e2 


This inequality states that the larger we take N the more probable is to find gN 
around (1/z) f fdv in a e neighborhood. The proof of this proposition use the decay 
of the auto-correlation function, for that reason the square integrable assumption on 
the observable /. 

Proof: First, notice that 


Ia^ 


gN 


m=l 

1 1 


N 1 - e-^/^ 


lA^ 
fdu 


<A^ 


(A.9) 


because v is cr-invariant. For the sake of clarity let us denote by / := fdu and let 
us define iF : N ^ R as follows. 


K{N) ■= 


1 


1 


N 1 - e-^/^ 
and note that for every 2 ; > 0 

lim KIN) = 

N^OO Z 

Observe that for any e > 0 we have 

/ 


9N{a) - 


>e> =iy 


gNiK)-KiN)f + KiN)f- 


f 


> e 


< V - K{N)f\ > e/2) + n {\K{N)f-f/z\ > e/2) . 

For the second summand in the right-hand side of the inequality, there is no 
probability. So, there exists a N*{e) such that for all N > N*, \K{N)f — f/z\ < e/2, 
and so that part eventually vanishes. While for the first summand, by the Chebyshev 
inequality we have that for all G N, 

!<,»(>) 


The variance oi gN can be rewritten in terms of the auto-correlation function, 

Var(5Ar) = J g%dv - gNdh^ = J glfdv - K‘^{N)p. 

We need to calculate / g%dh' and we obtain that 

/ -| OO OO n 

gjfdv = ifoa-)ifoP)d., 

m=lk=l d 
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1 OO OO 

Var(5iv) = 


^-mzlN^-kzlN 


m—1 k—1 

OO OO 


/ • / o (T™ ''dv - { I f du 


= e—- k). 

m—1 k—1 

where Cf stands for the auto-correlation function of the observable /, i.e., 


(A.IO) 


■= J f ■ f dn - fdiy^ . 


The expression for the variance of (/at given in Eq. (A.IO) can rewritten by changing 
the summation order, as follows 

- OO OO - OO 

Var((7Ar) = — 


1=1 n=l 


OO 

= KiN/2)-J2e-^^/^Cf{£) + -Q(0) 


n—1 

K{N/2) 


t=i 


Now, the assumption on the decay of correlations becomes cricial. Let us assume that 
the auto-correlation function for / behaves as Cf{£) = Co£~°‘ for any a > 0 and some 
Co > 0. In this case we can write the variance of pN as follows 

Var(g^) = K{N/2)^ 


1=1 


= K{N/2)^UUe-^/^) + . 

In the second equation we have rewritten the expression for the variance of gjv using 
the polylogarithm function (A.5). From the approximation (A.6) we are able to write 
down an asymptotic expression for the variance of for a sufficiently large N, 

Var( 5 ^) « K(iV/2)^r(l - a) + lc/(0)^^^. 

Finally, since K{N/2) ^ 1/z as A^ ^ oo, we easily see that Var((/Ar) 0 in the same 
limit, and specihcally 

Var(5w) «Cor(l-a)z“-2lV-“. 

The latter implies that we can chose an appropriate constant C such that for a 
sufficiently large N 


ffw(a) - 


f 


> e 


< 


C ■ r(l — a)z 


a-2 


Which concludes the proof the proposition. 


□ 
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Appendix B. The n-dimensional marginals for the expansion-modification 
system 


The 1-dimensional marginals for the E-M system can be understood as the relative 
frequency of occurrence of symbols (monomers) of every kind along a typical 
realization of the symbolic sequence (the polymer) generated by the E-M process 
at the stationary state. Eor the sake of simplicity here we will call “polymer” to an 
infinite symbolic sequence (i.e., a point in the space {0,1}^). 

The processes generating the infinitely long polymer of the E-M system, give us 
the way to obtain the marginals as follows. Eirst, let ap € {0,1}^ be an “initial 
polymer” (of length 0) and allow it evolve, coordinatewise, according to the E-M 
processes defined as 


J X with probability 1 — p 
( XX with probability p. 

These rules generate a sequence of “polymers” of finite size Hq i-?- ai i-> a 2 • • • that 
attain a stationary state for t —?■ oo. The E-M system is actually a discrete-time 
Markov process [21]. Let us call /o(t) and /i(t) the fraction of O’s and I’s in the 
symbolic chain a^ at time t G N. This means that in any subsequence of a^ of 
n symbols we will find typically nfo{t) zeroes and nfift) ones. We will obtain a 
recurrence relation for /q and fi. First notice that the chain a^+i is generated from 
a.t by means of a fraction of p expansions and 1 — p modifications. This means that 
a subsequence b of size n in at generate a sequence of size n(l + p) in at+i. The 
total symbols that are expanded in b is pn and the total symbols that are modified 
is (1 — p)n. From the total symbols that are “expanded”, a fraction /o of them are 
zeroes and a fraction /i of them are ones. Analogously from the total symbols that 
are “modified”, a fraction /o of them are zeroes and a fraction fi of them are ones. 
This means that the number of symbols that are zeroes in the new chain has two 
contributions: the number of O’s that were expanded 2np/o(t) plus the number of I’s 
that were modified n(l — p)fi{t) becoming zeroes in the new chain. An analogous 
statement follows immediately for fi{t) . Since the size of the new chain is (1 + p)n 
we have that the fractions fo{t) and fi{t) obey the following recurrence relations. 


s u . Wo(0 + (1 -p)/lW 

foit + ^) = -- 

+ 1) = + Wl(f) 


(B.l) 

(B.2) 


Notice that this defines a one-step Markov chain with probability matrix Mi 
{0,l}x{0,l}^[0,l], 


Ml 


1 

1+p 


2p 1-p \ 
1-p 2p ) ' 


(B.3) 


In the limit of t > oo the fractions /o and fi tend to the unique invariant probability 
vector for Mi. It can be seen that such a probability vector is given by 


Vi = 



(B.4) 


Next we can generalize the concept of fractions of zeroes and ones, to the relative 
frequency ft{^) with which a word w G {0,1}^ of size k occurs along a the polymer 
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Table Bl. Effective number of words of size two given by the expansion- 
modification process. This rules allows calculate the recurrence relations for the 
evolution of frequency of appearance of words of size two. 


Process 

Probability 


Initial words 


00 gives 

01 gives 

10 gives 

11 gives 

ee 


00 , 00 

00 , 01 

11 , 10 

11,11 

em 

p{i-p) 

00 , 01 

00 , 00 

11 , 11 

11 , 10 

me 

p{i-p) 

10 

11 

00 

01 

mm 

(i-pf 

11 

10 

01 

00 


SLt G {0,1}^ at time t. 

n—1 

ft{w) ■= lim ^I[w] (o-^(at)) (B.5) 

i=o 

where I is the indicating function, which takes the value 1 if a* belongs to the set [w] 
and take the value 0 if at does not belong to the set [w]. The set [w] called cylinder 
is defined as the set of all infinite sequences x G {0,1}^ with prefix w, i.e., 

[w] := |x G {0,1}^ : Xq = Wo,Xi =Wi,.. .,Xk-i = Wfe-i|. 

In other words, we calculate the frequency that a word of size k appears in a polymer 
by counting the number times such a word appears in a “window” of size k which 
slides along the polymer. Let us analyze the case k = 2. To obtain the stationary 
2-dimensional marginals, or equivalently, the stationary relative frequencies of words 
of size 2, we can obtain recurrence relations as in the 1-dimensional case. First notice 
that a word of size 2 generate words of sizes 2, 3 and 4 by means of expansions 
and modifications. We have four cases that should be considered. For example 
take the word a = 00. This word can be transformed under the E-M process into 
b = 0000,100,001,11. Each case corresponds to one of the four possibilities of the 
process: ee, em, me and mm (e stands for “expansion” and m stands for “modification”). 
First notice that in the case of two expansions (00 0000) we obtain, by sliding a 

window of size two, three times the word 00. However, this fact does not mean that 
pure expansions has the effect of triplicate the frequency of appearance of 00. For 
example, the block 00000 transforms under pure expansion into the block 0000000000. 
A direct counting shows, by sliding a window of size two, that the word 00 appears four 
times before the transformation and nine times after such a process. In general, the 
number of times that the word 00 appears in the block resulting from pure expansions 
of a block consisting of n zeros, is 2n-|-1. This means that the frequency of appearance 
of 00 is duplicated in the limit of an infinite polymer. In other words, the effective 
number of words that produces the word 00 under the process ee is two and are 00, 
00. A similar analysis shows that, in order to count correctly the number of times a 
word appears under the processes em, me, and mm we should take into account certain 
“rules” which are summarized in Table Bl. With these rules it is easy to see that the 
relative frequencies with which the words 00, 01, 10 and 11 appear after a process of 
expansions and modifications of an infinite polymer, are given by, 

/t+l(00) = ^ (^(2/ +pg)/t(00) + (/ + 2pq)M01)+pqMW) + q^MU)^ 
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/t+i(01) = ^(^pqft{00)+p^ft{01) + gVt(lO) + w/i(ll)) 

/t+i(10) = Y^(^M/t(00) + gVt(01)+pVt(10)+M/i(ll)^ 

/t+i(ll) = Y^(^'?Vt(00) +M/t(01) + {p"^ + 2pq)ft{10) + {2p‘^ +pq)ft{n)^. (B.6) 

The above recurrence relations can be recast into a more compact form by defining 
the probability vector p 2 (the time-dependent two-dimensional marginals) whose 
components are the relative frequencies at time t, 


P2{t) := (^ft+iiOO), ft+iiOl), ft+iilO), /i+i(ll)), 


in terms of which we rewrite Eq. (B.6) as a matrix equation, 

P2{t+ 1) =P2{t)M2 


where the matrix M 2 dehned as 


M 2 


l+p 


2p^ 

- 


+ pq 

V2pq 


pq 


pq pq q^ \ 
p^ q^ pq 

q2 p2 p2 _|_ 2pq 

pq pq 2p^ +pq J 


(B.7) 


In the limit of t —>■ 00, it has been proved in Ref [21] that the system reaches a 
stationary state. This means that the time-dependent probability vector p2(i) is well 
defined in the limit t —> 00. The probability vector at infinity is the left eigenvector 
of M 2 , which we denote by U2, associated to the (largest) eigenvalue A = 1, 


V 2 = lim P 2 {f), 

t—^OO 

and V 2 satisfies the equation 


V2 = V 2 M 2 , 

through which we calculated the 2-dimensional marginal given in Eq (26). 
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